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Abstract 

We consider a 1-parameter family of self-adjoint extensions of the 
Hamiltonian for a particle confined to a finite interval with perfectly re- 
flecting boundary conditions. In some cases, one obtains negative energy 
states which seems to violate the Heisenberg uncertainty relation. We 
use this as a motivation to derive a generalized uncertainty relation valid 
for an arbitrarily shaped quantum dot with general perfectly reflecting 
walls in d dimensions. In addition, a general uncertainty relation for 
non-Hermitean operators is derived and applied to the non-Hermitean 
momentum operator in a quantum dot. We also consider minimal un- 
certainty wave packets in this situation, and we prove that the spec- 
trum depends monotonically on the self-adjoint extension parameter. In 
addition, we construct the most general boundary conditions for semi- 
conductor heterostructures such as quantum dots, quantum wires, and 
quantum wells, which are characterized by a 4-parameter family of self- 
adjoint extensions. Finally, we consider perfectly reflecting boundary 
conditions for relativistic fermions confined to a finite volume or local- 
ized on a domain wall, which are characterized by a 1-parameter family 
of self-adjoint extensions in the (1 + l)-d and (2 + l)-d cases, and by a 
4-parameter family in the (3 + l)-d and (4 + l)-d cases. 
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1 Introduction 



The subtle differences between Hermiticity and self-adjointness of quantum me- 
chanical operators, which were first understood by von Neumann are rarely 
emphasized in quantum mechanics textbooks. This already affects the elementary 
textbook problem of a particle in a box jij. Almost exclusively, the students are 
taught to set the wave function to zero at the boundary, in order to ensure its con- 
tinuity. However, it is sufficient to guarantee that no probability leaks outside the 
box, i.e. that the probability current (but not necessarily the wave function itself) 
vanishes at the boundary. The resulting most general perfectly reflecting boundary 
condition contains a real- valued parameter that characterizes a family of self-adjoint 
extensions of the quantum mechanical Hamiltonian. In general, the wave function 
then does not go to zero at the boundary, and consequently the probability density 
to find the particle directly at the wall does not vanish. This is consistent with 
classical intuition of a ball bouncing off a perfectly reflecting wall. 

While the self-adjoint extensions of the Hamiltonian are certainly known to the 
experts, in the beginning of this paper we introduce them in a pedagogical manner, 
since they, unfortunately, seem not to constitute common knowledge in quantum 
mechanics. When re-doing the standard textbook problem of a particle in a box, we 
will find states of negative energy for a particle with only kinetic energy. Such states 
seem to violate the Heisenberg uncertainty relation. In order to clarify this issue we 
then derive a generalized uncertainty relation valid for a particle confined to an ar- 
bitrarily shaped region with general perfectly reflecting walls in d dimensions. This 
situation is relevant in the context of quantum dots. The most general perfectly 
reflecting boundary condition ensures that the component of the current normal to 
the reflecting surface must vanish. Again there is a 1-parameter family of self-adjoint 
extensions. The real-valued parameter that characterizes the boundary condition 
is a material-specific constant whose value could be determined experimentally for 
actual quantum dots. In particular, this parameter enters the generalized uncer- 
tainty relation. We also consider minimal uncertainty wave packets in a general 
quantum dot. Interestingly, when the self-adjoint extension parameter vanishes, a 
constant wave function of zero energy has Ap = and saturates the generalized 
uncertainty relation. While the purpose of our paper is to some extent pedagogical, 
to the best of our knowledge the generalized uncertainty relation for quantum dots 
has not been derived before. Perfectly reflecting boundary conditions for relativistic 
fermions have been investigated in the context of the MIT bag model [3.-;5|. Here we 
construct the most general boundary condition for relativistic Dirac fermions, which 
is characterized by a 1-parameter family of self-adjoint extensions in the (1 + l)-d 
case, and by a 4-parameter family in the (3 + l)-d case. Finally, we extend the 
discussion to domain wall fermions in (2 + l)-d and (4 + l)-d space-times. 

The paper is organized as follows. In section 2 we discuss a particle confined 
to a finite interval with general perfectly reflecting boundary conditions, and we 
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derive the corresponding generahzed uncertainty relation. In section 3 this relation 
is extended to an arbitrarily shaped quantum dot with perfectly reflecting walls 
in d dimensions. We also derive a general uncertainty relation for non-Hermitean 
operators and apply it to the non-Hermitean momentum operator. In addition, we 
construct the corresponding most general minimal uncertainty wave packet, and 
we prove that the spectrum varies monotonically with the self-adjoint extension pa- 
rameter. We also construct the most general boundary conditions for semiconductor 
heterostructures. In section 4, we extend the discussion to relativistic Dirac fermions 
in (1 -|- l)-d and (3 + l)-d, and in section 5 to domain wall fermions in (2 + l)-d and 
(4 + l)-d space-times. Finally, section 6 contains our conclusions. 



2 Particle in a 1-d Box with General Perfectly 
Reflecting Walls 

Let us consider a particle of mass m moving in the 1-d interval VL = [—L/2, L/2\. 
This problem has been discussed in the context of self-adjoint extensions in 
Other examples of quantum mechanical problems involving the theory of self-adjoint 
extensions are discussed, for example, in jof. We use natural units in which h = 1. 
For simplicity, we restrict the Hamiltonian to the kinetic energy operator 

H = I^ = -—dl (2.1) 
2m 2m ^ ^ 

The wave function \I'(x,t) gives rise to the probability current density 

t) = -^[^(x, tyd.^ix, t) - ty^ix, t)], (2.2) 

which together with the probability density p{x,t) = |\I'(x,t)p obeys the continuity 
equation 

dtp{x,t) + d,j{x,t) = 0. (2.3) 

In the following discussion, we can ignore the time- dependence of the wave function 
and simplify the notation to \E'(x). 



2.1 Spatial Boundary Conditions 

In order to guarantee probability conservation, one must demand that the probabil- 
ity current vanishes at the boundary, i.e. 

j{L/2)=j{-L/2) = 0. (2.4) 

The most general local boundary condition that implies this takes the form 

7(a;)^(a;) + d^^{x) =0, x = ±L/2. (2.5) 
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Indeed, one then obtains 



j{x) = ^[-^(x)*7(x)^(x) + 7(x)*^(a;)*^(x)] = ^ 7(0;) G R, x = ±L/2. 

(2.6) 

The two real-valued parameters 7(L/2) and 7(— L/2) characterize a 1-parameter 
family of self-adjoint extensions of the Hamiltonian at each of the two ends of the 
interval Q. In order not to break parity via the boundary conditions, we restrict 
ourselves to 

7(L/2) = -ii-L/2) = 7 G M, (2.7) 

such that 

7^(L/2) + a,^(L/2) = 0, -7^(-L/2) + 9,^(-L/2) =0. (2.8) 



2.2 Self-Adjointness of the Hamiltonian 

In order to investigate whether the Hamiltonian is indeed self-adjoint when the wave 
functions obey the boundary conditions eq. fl2.8p . we now consider 

1 /"^/^ 

(Xli^l^) = -7^ dx xixTdl^ix) 

2"^ J-L/2 

1 /"^/^ 1 

= ^ / dx d^xixTd^^ix) - — [x{xrd,^{x)]% 

= y dx dlxixY^ix) + — [dMxT'^ix) - x{xrd,^{x)]% 

= {^^J\H\xr + ^[^Mxr'^ix)-x{xr^,<^I{x)]%. (2.9) 

The Hamiltonian is Hermitean (or symmetric in mathematical parlance) if 

{X\H\^) = {H^xm = {Hxm = {m\xr, (2.10) 
which is indeed the case if 

[dMxT'^ix) - x{xrd,'^{xtJL/2 = 0- (2-11) 

The domain D{H) of the Hamiltonian contains the at least twice-differentiable 
square- integrable wave functions \E'(x) that obey the boundary condition eq. (12.81) . 
Using that condition, eq.f l2.1ip reduces to 

[d,x{L/2y + 7X(V2)*] - ^(- ^2) [9.x(-L/2)* - 7X(-V2)*] = 0. 

(2.12) 

Since \E'(L/2) and \1'(— L/2) can take arbitrary values, the Hamiltonian is Hermitean 
if 

lx{L/2) + d^x{L/2) = 0, -ix{-L/2) + d^x{-L/2) = 0, (2.13) 
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i.e. if the wave function x{^) c'-lso obeys the boundary condition eq. fl2.8p . Imposing 
this boundary condition also on x{^) imphes that the domain of H'^ coincides with 
the domain of H, D{H^) = D{H). Since H is indeed Hermitean when both 
and x{^) obey eq. fl2.8p . and since, in addition, D{W) = D{H), the Hamihonian is, 
in fact, self- adjoint. 

It should be noted that there is even a 4-parameter family of self-adjoint exten- 
sions of if Here we have encountered only two parameters, 7(L/2) and 7(— L/2), 
which we have reduced to one parameter 7 by demanding parity symmetry. The 
other two parameters of the 4-parameter family of self-adjoint extensions relate the 
values of the wave function and its derivative at x = L/2 to the corresponding values 
at a; = —L/2. Such a boundary condition violates locality and is thus not physically 
meaningful in the present context. 



2.3 Non-Hermiticity of the Momentum Operator 

In order to investigate whether the momentum operator p = —idx is self-adjoint or 
at least Hermitean, let us also consider 

rL/2 



J -L/2 
rL/2 

= I dx dxxixT'^ix) - I [x(x)*M/(a;)]% 

J-L/2 

= mP\xY-^[x{xy^{x)]%. (2.14) 

Hence, p would be Hermitean only if 

x(L/2)*M/(L/2) = x(-L/2)*^(-L/2). (2.15) 

There is no reason why this should be the case when '^{x) and x(x) obey the 
boundary condition eg. (12.81) . Hence, in the domain D{H) of the Hamiltonian, the 
momentum operator p is not even Hermitean, and thus certainly not self-adjoint. 

The only self-adjoint extension of the momentum operator on a finite interval 
results from the boundary condition 

^(L/2) = A^(-L/2), A G C. (2.16) 

Inserting this relation in eq. fl2.15p . we obtain 

x{L/2yX^{-L/2) = x(-L/2)*xI/(-L/2) ^ x{L/2) = ^x(-L/2). (2.17) 

If \I'(x) obeys eq. fl2.16p and x(x) obeys eq. fl2.17p the operator p is Hermitean (i.e. 
symmetric). The domain D{p) contains those at least once-differentiable square- 
integrable wave functions \I'(x) that obey eq. fl2.16p . The domain D[p^), on the 
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other hand, contains the corresponding wave functions x{^) that obey eg. (12. 171) . 
The operator p is self-adjoint only if D{p) = D[p^) which implies 

A = ^ A = exp{i9) 

= exp{z9)^{-L/2), x{L/2) = exp{z9)x{-L/2). (2.18) 

Hence, the momentum operator is self-adjoint only if the probability density is a 
periodic function, i.e. piL/2) = |^(i:/2)|2 = \^{-L/2)\^ = p{-L/2). Since non- 
local periodic boundary conditions make no physical sense in the present context, 
and since the wave functions in the domain of the self-adjoint Hamiltonian obey 
the boundary condition eq. (l2.8p . but not eq. (l2.18p . in this case the momentum 
operator is neither Hermitean nor self-adjoint. Hence, we must conclude that, in the 
present context, momentum is not a physical observable. This indeed makes sense 
for a particle that is confined to a finite region of space. After all, a momentum 
measurement would put the particle in a momentum eigenstate, which also exists 
outside the box and would therefore require infinite energy. An alternative point 
of view is taken in js], [ij where the infinite potential in the energetically forbidden 
region is approached as a limit of a large but finite potential. 



2.4 Energy Spectrum and Energy Eigenstates 

Let us now consider the energy eigenstates and the corresponding energy eigenvalues 
for the particle in a box with general reflecting boundary conditions parameterized 
by 7. First, we consider positive energy states of even parity, i.e. 

^n{x) = Acos{knx), ^ = n = 0, 2, 4, . . . , OO. (2.19) 

2m 

The boundary condition eq. (12.81) then implies 

7cos(fc„L/2) - A;„sin(fc„L/2) = ^ = tan(A;„L/2). (2.20) 

Interestingly, for 7 = there is a zero-energy solution with = and a constant 
wave function '^^{x) = ^/l/L. 

Similarly, the positive energy states of odd parity take the form 

^n{x) = Asm{knx), En = ^, n = 1, 3, 5, . . . , OO, (2.21) 

and must obey 

7sin(A;„L/2) + A;„ cos(A;„L/2) = ^ = - cot(A;„L/2). (2.22) 
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In this case, a zero-energy solution exists for 7 = — 2/L with the wave function 
= ^/V2/I? X. While this solution emerges from ^i(x) = Asm.{kix) in the 
limit 7 — > — 2/L, it also follows directly from the zcro-cncrgy Schrodinger equation 
dl-^{x) = 0, and the boundary condition (-2/L)*(L/2) + d^^{L/2) = 0. 

Next, let us consider eigenstates of negative energy and even parity. In that case, 
the wave function takes the form 

^^[x) = Acosh(xa;), Eo = - — , (2.23) 

2m 

and must obey 

7COsh(xL/2) + xsinh(xL/2) = =^ - = -tanh(xL/2). (2.24) 

X 



Again, for 7 = one recovers the zero-energy state '^o{x) — ^/l/L. 

Finally, we consider the negative energy eigenstates with odd parity, i.e. 



x^ 



^i{x) = ^sinh(xx). El = - — , (2.25) 



which must obey 



7sinh(xL/2) + xcosh(xL/2) = ^ - = -coth(xL/2). (2.26) 

X 

In this case, for 7 = —2/L we recover the zero-energy eigenstate ^i(x) = ^JV2JU' x. 

The spectrum as a function of 7 as well as the corresponding wave functions 
are illustrated in figure 1. For 7 = 00, we recover the standard textbook case with 
'^{L/2) — ^(— L/2) = 0. In that case, the energy spectrum takes the famihar form 

£;„(7 ^ 00) = ^^^^i^^ , n = 0,l,2,...,oo. (2.27) 
As 7 decreases to zero, the energy eigenvalues decrease such that 

K(7 = 0) = ^^, n = 0,l,2,...,oo. (2.28) 

For —2/L < 7 < there is one negative energy state, and for 7 < —2/L there are 
even two negative energy states, which reach negative infinite energies in the limit 
7 —00, i.e. 

£^0,1(7 ^ -00) ^ ^ -00. (2.29) 

In the limit 7 —00, the wave functions of the negative energy states reduce to 
(5-function-like structures localized at the boundaries. The rest of the spectrum is 
exactly as in the standard textbook case (i.e. 7 = 00), namely 

K(7^-oo)^ , n = 2,3,4,...,oo. (2.30) 
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arctan(yL/2) 




Figure 1: Top panel: Energy spectrum of a particle in a 1-d box as a function of the 
self-adjoint extension parameter 7. The x-value represents arctan(7L/2) G [— f , f ]; 
which corresponds to •j E [— cxd,cxd]. The y-value represents the energies En (with 
n = 0, 1,2,3, 4j in units 0/ 7r^/(2mL^). Bottom panel: The wave functions "^nix), 
X G [—■§,■§], (with n = 0,l,2,3j, for 7 = — 00, — |;, 0, 00. The sharp peaks in 
the n = and n = 1 states at '-f = —00 represent 6 -function-type wave functions 
of negative infinite energy localized at the boundaries. Except for these states, the 
energies and wave functions at 'y = —00 are the same as those 7 = 00. 
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The existence of negative energy states is somewhat counter-intuitive. In partic- 
ular, since the boundary conditions are perfectly reflecting, one may expect that the 
particle cannot bind to the wall. However, quantum mechanics does indeed allow 
the existence of walls which are perfectly reflecting for positive energy states and, 
at the same time, "sticky" for negative energy states. Since negative energy states 
have E = {p^ /2m) < 0, they seem to violate the Heisenberg uncertainty relation 
AxAp > i. Obviously, any state confined to the box has Ax < L/2. Hence, the 
Heisenberg uncertainty relation seems to suggest that Ap > 1/L, which would be 
in conflict with (p^) < 0. To resolve this puzzle, one must realize that the Heisen- 
berg uncertainty relation was derived for an infinite volume, and thus needs to be 
reconsidered in a finite box with perfectly reflecting boundary conditions. 

2.5 Uncertainty Relation for a 1-d Box with Perfectly Re- 
flecting Boundary Conditions 

Let us derive a variant of the Heisenberg uncertainty relation, taking into account 
the boundary condition eq. fl2.8l) . We follow the standard procedure by constructing 
the non-negative integral 

rL/2 

1=1 dx\d^^{x)+ax^{x)+ (3^{x)\^ >Q. (2.31) 

J -L/2 

Here a G M and (5 = fir + iPi C are parameters to be varied in order to minimize 
/ and thus derive the most stringent inequality. Using partial integration as well as 
the boundary condition eq. fl2.8l) . it is straightforward to obtain 

/ = (p2) _ + ^ p-i + ct2^a;2^ + 2a/3^(x) + 2/3ip + aa-h + /3rC, (2.32) 
where we have introduced 

a = ^[p{L/2) + pi-L/2)l 
h = ^[p{L/2 + p{-L/2)l 

c = p(L/2)-p(-L/2), p{±L/2) = \^>{±L/2)\\ (2.33) 
We have also defined 

rL/2 

p=^ dx ^{xy{-id^)^{x). (2.34) 

J -L/2 

If p = —idx were a self-adjoint operator, p would simply be the momentum ex- 
pectation value. However, since p is not self-adjoint in this case, the corresponding 
expectation value is in general complex. In particular, p is not the expectation value 
of any observable associated with a self-adjoint operator. 
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By varying a, /3r, and /3j in order to minimize /, one obtains 



dl 

— = -1 + 2a{x^) + 2l3r{x) + a = 0, 
oa 

dl c 

— =2f3r + 2a{x) + c = /3, = -a{x) - -, 

dl 

— = 2A + 2p = ^ A = -p, (2.35) 



which imphes 



a = i±^^, {AxY = {x')-{x)' = {{x-{x)Y). (2.36) 

Inserting these results for a, /3r, and back into the expression for /, eq. fl2.32p . one 
obtains 

I={p')-f-{ -&-^>0, (2.37) 



2Ax J 4 
which imphes the generahzed uncertainty relation 

<P>^^+( 2Ax ) +"+4- 

In the infinite volume limit, the wave function vanishes at infinity, such that a = 
b = c = 0. Furthermore, the momentum operator would then be self-adjoint with 
(p) = P- Hence, in the infinite volume limit, the inequality fl2.38p reduces to the 
standard Heisenberg uncertainty relation 

{p') > {pr+ (^X ^xAp > i, {Apr = {p')-{py = {{p-{p)r). (2.39) 



Interestingly, thanks to time-reversal invariance, the eigenfunctions of the par- 
ticle in the 1-d box are real-valued and thus p = 0. Since H = p^/2m, for energy 
eigenstates the generalized uncertainty relation takes the form 

2mK = ip') > (^^1^^) ' + (2.40) 

The zero-energy eigenstate \E'o(a;) = a/V^ for 7 = has a = l, 6 = c = 0, such 
that the uncertainty relation reduces to 2mEo > 0, which is indeed satisfied as an 
equality. This means that this state represents a minimal uncertainty wave packet. 
The zero-energy state \E'i(x) = ^JvIJU" x for 7 = — 2/L, on the other hand, has 
a = 3, 6 = — 12/L^, c = 0, as well as Ax = a/3/20 L, such that the inequality 
reduces to 2mEi > — 16/(3L^), which is again satisfied, but not saturated as an 
equality. Hence, in this case, the state does not represent a minimal uncertainty 
wave packet. 
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2.6 Constructing a Wall with 7 < oo 



It is interesting to ask whether perfectly reflecting walls with 7 < 00 are just a math- 
ematical curiosity or whether they can also be constructed physically. Of course, 
it is clear that, independent of the value of 7, a perfectly reflecting wall is always 
a mathematical idealization. Any real wall will eventually be penetrable if it is hit 
by a sufficiently energeticparticle. In this context it may be interesting to men- 
tion the MIT bag model [sHsl, in which quarks are confined by restricting them to 
a finite region of space with perfectly refiecting boundary conditions. While this 
model is quite successful in modeling the confinement of quarks, in the fundamen- 
tal QCD theory confinement does not arise through perfectly refiecting walls but 
through confining strings. When a dynamical QCD string connects a very energetic 
quark-anti-quark pair, the string can break by the dynamical creation of further 
quark-anti-quark pairs. Since perfectly refiecting walls do not even exist in the con- 
fining theory of the strong force, they certainly do not exist in condensed matter 
either. However, when low-energy particles hit a very high energy barrier, it acts 
effectively like a reflecting wall. In this sense, the self-adjoint extension parameter 7 
can be viewed as a low-energy parameter that characterizes the reflection properties 
of a very high energy barrier. 

Since 7 appears as a natural mathematical parameter characterizing a perfectly 
reflecting wall, there is no reason to expect that it cannot be physically realized. 
To show this explicitly, we now construct a wall with an arbitrary value of 7 as a 
limit of square-well potentials^ Let us consider the potential illustrated in flgure 
2, which consists of a perfectly reflecting wall at a; = with the standard textbook 
value 7 = cxD, and a very narrow and very deep square- well potential of size e > 
and depth —Vq < next to it. The wave function then takes the form 

^(x) = Asin(gx), < x < e, E = - Vq. (2.41) 

2m 

For X > e, the wave function is determined by enforcing continuity of the wave 
function itself and its derivative at x = e. When e — i- 0, we can determine the 
resulting value of 7 from 

- 7*(e) + (9^^(e) = 7 = gcot(ge). (2.42) 

In order to keep 7 flxed as e — )• 0, we must hence let q go to inflnity such that 

TT 2 , , 

q=- 7- 2.43 

2e TT 

For all states of flnite energy E, this is achieved by sending Vq = /2m — )■ cxd, in 
such a way that 

^We like to thank F. Niedermayer for suggesting this construction. 
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Figure 2: A deep and narrow square-well potential V{x) with depth — Vq — > — oo and 
width e — mimics a boundary with non-standard self-adjoint extension parameter 
7 < oo. The wave function and its first derivative are continuous at x = e. In 
the example shown here, 7 < 0, such that a bound state is localized at the boundary. 

It goes without saying that a perfectly reflecting wall with an arbitrary value of 7 
can also be constructed in many other ways, for example, by using an attractive 
(5-function potential next to the wall. 

2.7 Experimental Determination of 7 

It is natural to ask how one can determine the value of 7 for some perfectly reflecting 
wall that can be investigated experimentally. For example, for a planar homogeneous 
perfectly reflecting wall, the material-specific parameter 7 can be determined from 
the scattering phase shift 6{k) of an incident plane wave that propagates perpen- 
dicular to the surface. Assuming that the incident wave propagates in the negative 
x-direction in the region a; > 0, and scatters off a perfectly reflecting wall at a; = 0, 
the wave function takes the form 

"^{x) = exp{—ikx) + Rexp{ikx), x > 0. (2-45) 

Imposing the boundary condition 

-7^(0) + 9^^(0) = 0, (2.46) 
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one obtains 

- 7(1 + R)- ik(l -R) = =^ R = exp(i6(k)) = _I±1^. (2.47) 

7 — i/c 

Hence, by measuring the phase shift 

S{k) = 2 arctan(A;/7) + tt, (2.48) 
one can determine the material-specific self-adjoint extension parameter 7. 

3 Uncertainty Relation for a Quantum Dot with 
General Perfectly Reflecting Boundary Condi- 
tions 

In this section we consider an arbitrarily shaped rf-dimensional region Q with general 
perfectly refiecting walls at the boundary dQ. This may be viewed as a model for 
a quantum dot, in which electrons are confined inside a finite region of space. As 
illustrated in figure 3, Q may have multiple disconnected boundaries. For a math- 
ematical exposition of boundary value problems for operator differential equations 
we refer to [lOj]. Let us consider the Hamiltonian 

H = ^ + Vix) = --^A + V{x), (3.1) 
2m 2m 

where V{x) is a non-singular potential. In d dimensions, the continuity equation for 
local probability conservation takes the form 

dtp{x,t) + V ■j{x,t) = 0, (3.2) 
with the probability density and current density given by 

p{x,t) = \^{x,t)\\ j{x,t) = ^ [*(f,t)*W(f,t) - W(f,t)*^(f,t)] . (3.3) 

Again, in the following discussion the time-dependence is not essential and we sim- 
plify the notation to \E'(a;). 



3.1 Spatial Boundary Conditions 

As in the 1-d case, we demand that no probability may leak outside the region Q. 
This is ensured by requiring that the component of the probability current density 
normal to the surface vanishes, i.e. 

n{x) ■ j{x) = 0, xe dfl, (3.4) 
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Figure 3: A region Q of a quantum dot with two disconnected boundaries dQ. The 
shaded region is energetically forbidden. The vector n is normal to the boundary. 

where n{x) is the unit-vector normal to the surface dVL at the point x. The most 
general local boundary condition that ensures this takes the form 

7(f)^'(f) + n{x) ■ V^(f) = 0, X e dn, (3.5) 

and one then indeed obtains 

n(f) ■ = ^[-^(f)*7(x)^(x,t) +7(f)*^(x,t)*^(f,t)] = 0, xedn, (3.6) 

which again implies 7(0;) G M. In general 7(x) will depend on the position x G dVL 
on the boundary. In a real system such as a quantum dot, 7(0;) is a material-specific 
parameter, to be determined experimentally. 
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3.2 Self-Adjointness of the Hamiltonian 



Let us again convince ourselves that the Hamiltonian endowed with the boundary 
condition eq. fl3.5p is indeed self-adjoint. For this purpose, we consider 



{X\H\^) 



2m 



A + V{x) 



2m 



an 



d^x 



2m 

dn ■ x(f)*V\l'(x) 
1 



+ 



1 
2m 



x{xy\^{x 



dn 



Vx(x)*^(^) - x(f)*V^(f) 



do. 



mH\xY + 7^ dn- VxixY^ix) - x(x)*V^(x) 
Jan 



Hence, the Hamiltonian is Hermitean if 



dn ■ 



an 



Using the boundary condition eq.f l3.5p . the integral in eq. (l3.8p reduces to 



d'^-^x 



an 



n{x) ■ Vxi^y + 7(x)x(a;)* \l'(x) = 0. 



(3.7) 



(3.8) 



(3.9) 



Since \E'(x) itself can take arbitrary values at the boundary, the Hermiticity of H 
requires that 

n{x) ■ Vx(f) + 7(f)*x(^) = 0. (3.10) 

For 7(0;) G M, this is again the boundary condition of eq. fl3.5p . which ensures that 
D{H^) = D{H), such that H is indeed self-adjoint. In complete analogy to the 
1-d case, it is easy to convince oneself that the momentum operator p = —iV is 
not even Hermitean in the domain D{H) of the Hamiltonian, which contains those 
twice-differentiable and square-integrable wave functions that obey the boundary 
condition eq. fl3.5p . 



3.3 Generalized Uncertainty Relation 



Mathematical investigations of generalized uncertainty relations can be found, for 
example, in UJ, ll2|. In this subsection, we derive a generalized uncertainty relation 
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for an arbitrarily shaped quantum dot in d dimensions, by considering the non- 
negative integral 

d'^x |V^(x) + ax^{x) + > 0. (3.11) 



n 

Again, a G M and (5 = P^ + iPi & C'^ are parameters to be varied in order to minimize 
/. Using the boundary condition eq.f l3.5p after performing some partial integrations, 
one finds 

I = {^^)-ad + (3^ + (3f + a^{x^)+2a/3r-{x)+2(3i-f+a{n-x)-{-f)+(3r-{n). (3.12) 
Here we have defined 

{n ■ x) = dn ■ xp(x), 

JdQ 

(7) = / d'-'x i{x)p{x), 
Jan 



(n) = / dnp{x), p{x) = \^{x)\\ (3.13) 
Jdn 

In analogy to the 1-d case, we have also introduced 

j d'^x ^(f)*(-iV)^(f). (3.14) 

If p = — iV had been a self-adjoint operator, p would be the momentum expectation 
value. However, in this case, p is again not the expectation value of any observable 
physical quantity. 

By varying a, and one obtains 



da 



d + 2a(x ^) + 2/3,. ■ (f) + (ra ■ f ) = 0, 



dl ■ - in) 

— = 2Pi+ 2a{x^) + {n^) = 0^/3,, = -a{x) - 
dpi 2 

dl — - - 

— = 2Pj + 2pi = ^ A = -p, (3.15) 

8131 

which then implies 

« = ^^^^^^^^^y^^, {^xf = {x')-{xf = {{x-{x)f). (3.16) 

— * — # 

Again, by inserting the results for a, /S^, and Pi back into eq. fl3.12p . one obtains 

/ = {p=> -r - (i±M_<^y _ (-,) _ <|! > 0, (3.17) 
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which finally implies the generalized uncertainty relation 

{p-> > f + ( '^^"^ 'it" ' + (7) + ^- (3.18) 

It is not obvious that eq.f l3.18p is invariant under a spatial translation x — )■ x' = 
x + d. While (p^), p, Ax, (7), and (n) are translation invariant by construction, (x) 
and (n ■ x) are not. It is thus reassuring to realize that the combination 

(n) ■ (x') — (n ■ x') = {n) ■ ((x) + d) — {n ■ x) — {n ■ d) = {n) ■ (x) — (n ■ x), (3.19) 

that enters the generalized uncertainty relation is indeed translation invariant. 

All quantities entering the uncertainty relation, except p, are directly related to 
physically observable quantities. However, since the momentum is not a self-adjoint 
operator and hence not observable in a finite volume, despite the fact that it is 
mathematically completely well-defined, the quantity p seems not to be physically 
measurable. In that case, the question arises how the generalized uncertainty rela- 
tion should be interpreted physically. Since momentum cannot even be measured 
inside a quantum dot, (Ap)^ = (p ^) — p can obviously not be interpreted as the 
uncertainty of momentum. Still, (p^) determines the energy of an eigenstate of the 
free particle Hamiltonian with V{x) = 0. Since time-reversal invariance guarantees 
that energy eigenstates have real-valued wave functions, for those states we know 
that p = 0. In that case, the generalized uncertainty relation reduces to 

^ ,^9, f d + (n) ■ (x) — (n ■ x)\'^ , , (n)^ 
2mE„ = if) > ^ ^ ^ ^-j + (7) + (3.20) 

which indeed contains measurable physical quantities only. 

3.4 General Uncertainty Relation for non- Hermit ean Oper- 
ators 

For two general self-adjoint operators A and B, the uncertainty relation is sometimes 
quoted as AAAB > ^\{[A, B])\. However, a more stringent form of the inequality 
is given by 

AAAB >\{AB)\ = \{AB) - {A){B)\, A = A-{A), B = B-{B). (3.21) 

In the case discussed before, the momentum operator is not even Hermitean. Hence, 
the question arises, how the uncertainty relation looks like for a general pair of non- 
Hermitean operators A ^ A^ and B 7^ B"^. In that case, it is natural to define 

{AAf = {A^\A^) = (A^A) = (A^A) - {A^){A) = (A^A) - > 0. (3.22) 
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Here we have used (A^) = {A)*. Using the Schwarz inequahty, we then obtain 

> \{A<I/\m)\^ = \{A^B)\^ = \{A^B) - {A^){B)\^, (3.23) 

such that 

AAAB > I {A^B) - (A^) (B) \ . (3.24) 

Let us compare the general uncertainty relation fl3.24p for non-Hermitean opera- 
tors with the generalized uncertainty relation f l3.18p that we obtained for a particle 
confined to the region Q. In that case, we can identify A with the self-adjoint op- 
erator X and B with the non-Hermitean operator p = —iV. The uncertainty of the 
momentum is then given by 

{ApY = {p^.p)-\{p)\'. (3.25) 
Since p^ ^ p, in this case, p^ ■ p ^ p"^. In particular, we obtain 

{p^ -p) = (p^Ip^) = / d'^x V^(x)* ■ V^(f) 

Jn 

= - d'^x ^(x)*A^(f) + dn- *(x)*V^(f) 

Jo. J do. 

2\ / 



= {p') - / d^-^x 7(f)^(f)*^(f) = {p') - (7). (3.26) 
J an 

Furthermore, decomposing (p) into real and imaginary parts, we find 

(p) = 3? I d'^x ^(f)*(-2V)^(f) +i3 f d'^x ^(f)*(-iV)^(f) Un). (3.27) 

Here we have used 

$5 [ rf'^x ^(f)*(-iV)^(f) = — ! d'^x\^{xY{-iV)^{x) -^{x){iV)-^{x) 
Jn 2« L 

= -\ I dn^{xr^{x) = -Un), (3.28) 
^ Jan ^ 

and hence we obtain 

(Ap)2 = (j72)-(7)-/-^. (3.29) 
The generalized uncertainty relation f l3.18p can thus be rewritten as 

AxAp > (n) ■ (x) - (3.30) 

On the other hand, in this case the general uncertainty relation for non-Hermitean 
operators (I3.24p takes the form 

AxAp > |(x -p) - (x) ■ (3.31) 
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Introducing x ■ p = ^{x ■ p) and using 

■ (-zV)^(f) - ^{x)x ■ (zV)^(x)* 

= - I rf'^x \[^(f)*V ■ - - / ciri ■ 

2 in 2 

= i(rf-(n-x)), (3.32) 

one obtains 

{x ■ p) — (x) ■ (p) = X ■ p — (x) ■ p + -{d — {n ■ x) + (n) ■ (x)), (3.33) 
such that 



AxAp > \j {x ■ p — {x) ■ pY + -{d — {n ■ x) + (n) • (x))^. (3.34) 

Hence, unless x ■ p = (x) ■ p, the general uncertainty relation for non-Hermitean 
operators (13.241) is more stringent than the generalized uncertainty relation (I3.18p . 
Since energy eigenstates have a real-valued wave function, for them x ■ p = and 
p = 0, such that both inequalities are then equivalent. 



3(a; ■ p) 



i 

27. 



d'^x 



3.5 Minimal Uncertainty Wave Packets 

It is interesting to ask which wave functions saturate the generalized uncertainty 
relation fl3.18p . and thus satisfy it as an equality. First of all, when 7(x) = 
everywhere at the boundary, the constant wave function \I'(x) = 1/a/K, where V 
is the volume of fi, is an energy eigenstate of zero energy. In that case, we obtain 
p = 0, (7) = 0, as well as 

(n ■ x) = — I dn ■ X = — I d'^x V ■ x = d, 
(n) = 1 / dn = ^ d'^x VI = 0, (3.35) 



V Jan V 

such that the inequality (13.181) then reduces to 2mE = (p^) > 0. Hence, the zero- 
energy state indeed saturates the inequality. In this sense, it can be viewed as a 
state of minimal uncertainty. Of course, one should not forget that, since in this 
case momentum is not even a physical observable, (p^) cannot be interpreted as the 
uncertainty of a momentum measurement. 

Are there other minimal uncertainty wave packets beyond the zero-energy state 
with a constant wave function that exists for 7(x) = 0? It is clear by construction. 
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that the inequahty f l3.18p can be saturated only if the integrand in eq. (13. lip vanishes, 
i.e. if 

V^(f) + af^(f) + = 0. (3.36) 

Using the boundary condition eq.f l3.5p . for points x G dQ this imphes 

-7(x) + an{x) ■ x + n{x) ■ f3^ \l/(x) = ^ 7(x) = n{x) ■ {ax + (3). (3.37) 

Since 7(x) G M, we must further demand /3j = —p = 0. This imphes that 7(0;) must 
have a very pecuhar form at the boundary, which will generically not be the case. 
Thus, in general, there are no minimal uncertainty wave packets in a finite volume. 
In the bulk eq. fl3.36p is satisfied just by the standard Gaussian wave packet 

^(f) = Aexp , (3.38) 

which is known to saturate the Heisenberg uncertainty relation in the infinite volume. 
We then also have 

d+{n) ■ (f ) -{n-x) - (n) 
" 2(A^)2 ' ^ " ~ ~Y' ^ ^ ^ 

such that 7(x) must satisfy 

/ (n)\ (fi)^ 

7(f)=n(f)-U(f-(f))- Yj (^) = a((n-f)-(n)-(f))-^. (3.40) 

Only if 7(x) happens to be such that the Gaussian wave packet automatically sat- 
isfies the boundary condition, it remains a minimal uncertainty wave packet in the 
finite volume. 



3.6 7-Dependence of the Energy Spectrum 

In this subsection we assume that 7(x) = 7 G M is a constant independent of the 
position X G dQ on the boundary. We then ask how the energy spectrum changes 
with 7. A similar calculation for the Dirac operator in a relativistic field theory was 
performed in [l3[. Let us introduce the energy eigenvalues En and the corresponding 
wave functions ^^(a^), i-e. 

H^nix) = (-^A + V{x)^ ^„(f) = En^nix). (3.41) 

Both En and "^nix) depend on 7 via the boundary condition 

7^„(x) + n{x) ■ V^„(f) = 0, f G dQ. (3.42) 



20 



The 7-dependence of the energy spectrum follows from 

Jn 

d'^X [d^^n{xyH^n{x) + H^n{xyd^^n{x)] 



+ — / dn- V^n{xyd^^n{x)-^n{xYVd^^n{x) 
2"^ Jan L 

= En [ d'^x 
Jn 

+ d'-'x[-^^n{xyd^^n{x) + ^n{xyd^{^^n{x))] 

Jan 

= / d'^-'x ^„(f)*^„(f) = ^ [ d^-'x pn{x) > 0, (3.43) 

Jon Jan 

which shows that the spectrum is monotonically rising with 7. 

Indeed, in the 1-d case, the spectrum illustrated in figure 1 is monotonic in 7. 
In that case, eq.f l3.43p reduces to 



d.,Er, = — [Pn{L/2) + Pni-L/2)] . 



2m 



(3.44) 



Let us explicitly verify this equation for the positive energy eigenstates of even parity 



7 

\E'„(x) = y4cos(fc„x), En = -^, — = tau 

2m kr, 



In this case, one obtains 



d^E. 



dry kfi 2 k'fi cos^(A;„L/2) 



m 



m[knL + sin(/c„L)] 



as well as 



p„(L/2) + p„(-L/2) = 2\A\^ cos^ 



The normalization condition for the wave function implies 

1^1 



knL 1 
^ dx cos^ = I^Ptttt (^"-^ + sin(A;„L)) = 1, 



such that indeed 



1 2k 



2 ^nL „ 
COS — — = O^En. 



(3.45) 

(3.46) 
(3.47) 

(3.48) 
(3.49) 



It is straightforward to repeat this check for states of odd parity or negative energy. 
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3.7 General Boundary Conditions for Heterostructures 



Semiconductor heterostructures such as quantum dots, quantum wires, and quan- 



tum wells are separated into regions with different effective electron masses [14 
Until now we have considered quantum dots consisting of a single region isolated 
from an energetically forbidden environment. In this subsection, we construct the 
most general condition at the boundary dfl separating two regions, fli and fin as 
illustrated in figure 4, with different effective electron masses mi and mn- In the 
two regions the Hamiltonian is then given by 




Figure 4: Two regions Qi and fin in o quantum hetero structure are separated by a 
boundary dQ with the unit-vector n normal to the surface. An electron may have 
different effective masses mi and mn in the two regions. 
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For points at the boundary dQ, the Hermiticity condition takes the form 



nix] 



2mi 
1 

2mn 



(3.51) 



The self-adjointness condition can now be expressed as 



r(f) 



1 V r(x) e GL(2,C), (3.52) 



which turns the Hermiticity condition eg. (13. 51 p into 



r2i(^)Xi(^)* - rii(x)-^n(f) ■ Vxiix)* + -^n{x) ■ Vxn(^)* 

zmi Zniji 



*n(x) + 



r22(^)Xi(^)* - ^i2{x)- — n{x) ■ Vxiix)* - Xnix)* 

Znii 



In order to satisfy this relation we must thus demand 



2mn 



Xii(^) 
^n{x) ■ Vxn(x) 



r22(^) 



TAxY 



T2i{xy rn(x)* ) \ ^n(f).Vxi(f 



n{x) ■ V^'n(f) = 0. 

(3.53) 



XI (^) 



(3.54) 

Self-adjointness requires equality of the domains, D{W) = D{H), which thus implies 



1 



r22(x) -ri2(f) 

rii(f)r22(x) - ri2(f)r2i(f) V -^2i{x) T^^{x) 
r22(x)* -ri2(x)* 

-ToAxY Tu(xY 



(3.55) 



This condition is satisfied by a 4-parameter family of self-adjoint extensions which 
satisfy 

Tij{x) e exp{t9{x)) R, Tuix)T22ix) - ri2(x)r2i(x) = exp(22^(f)), (3.56) 
i.e. r(x) is real with determinant 1, up to an overall phase exp(i^(a;)). 

Using this form of r(a;), it is straightforward to show that self-adjointness implies 
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probability current conservation, i.e. 



in{x) ■ ji{x) = 



— 


n{x) ■ 


= 


n{x) ■ 


+ 


n{x) ■ 


+ 


\^e{x 


+ 


■ 




n{x) ■ 




inix) 



1 




2mi 




1 








1 


V 


2mn 





1 r 

— ^'i(x)*V^i(f) - V^i(x)*^i(x) 
2mi L 

L(^)T22(x)-r2l(f)Ti2(f)] 

,(f)T2i(f)-r22(f)Tn(f)] 
)|2[rn(f)T2i(f)-r2i(f)Tn(5)] 
-^V^n(x)|' [ri2(x)T22(x) - r22(x)Ti2(x)] 

[^n(x)*V*n(x) - V^n(x)*^n(x) 
2mn L 

Jn(x). 



= in{x) ■ j-a{x). (3.57) 

Here we have used again that r(a;) is real with determinant 1 up to the overall phase 
exp(z6'(a;)), which implies 



rii(f)T22(x) - r2l(x)Ti2(f) = 1, 

ri2(x)T2l(f) - r22(x)Tii(f) = -1, 

rii(f)T2i(f) - r2i(f)Tn(x) = 0, 

ri2(x)T22(x) - r22(x)Ti2(f) = 0. 



(3.5J 



4 Reflecting WaUs for Relativistic Fermions 



In this section we consider perfectly reflecting boundary conditions for relativis- 
tic fermions. Such boundary conditions have been introduced in the MIT bag 
model to mimic the confinement of quarks and gluons inside hadrons. Here we 
construct the most general perfectly reflecting boundary conditions for relativis- 
tic Dirac fermions both in one and in three spatial dimensions, as well as for the 
resulting non-relativistic fermions described by the Pauli equation. 



4.1 Reflecting Boundary Conditions for 1-d Dirac Fermions 



General boundary conditions for a Dirac particle in a box have been investigated in 



15[. Here we consider a Dirac fermion moving on the positive x-axis with a perfectly 
reflecting wall at a; = 0. The corresponding free particle Hamiltonian is then given 

by 

H = apc + (5mc\ a=(^J o)' ^=(o -l)" ^^"^^ 
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The Hamiltonian acts on a 2-component spinor \l/(a;,t), and the corresponding con- 
tinuity equation 

dtp{x,t) + dj{x,t) = 0, (4.2) 

is satisfied by 

p(x, t) = *(a;, t)"^^(a;, t), j{x, t) = c^(x, t). (4.3) 
Let us investigate the Hermiticity of the Hamihonian on the positive x-axis 

POO 

Jo 

= dx { [-caid^ + (3mc^] x{x)]^ ^{x) - icx(0)"^a^(0) 
Jo 

= {%H\xY-icx{Q)^a^{Q), (4.4) 
which thus leads to the Hermiticity condition 

X(0)"^a^(0) = 0. (4.5) 
We now introduce the self-adjoint extension condition 

^2(0) = A^i(O), AeC, (4.6) 

which reduces eq. fl4.5p to 

X(0)^f? i V(0) = [Xi(0)*A + X2(0)*]vl>i(0) = ^ X2(0) = -A*xi(0). (4.7) 



In order to guarantee self-adjointness of if, i.e. D{H) = D{W), we must request 

A = -A*, (4.8) 

i.e. A must be purely imaginary. For 1-d Dirac fermions, there is thus a 1-parameter 
family of self-adjoint extensions that characterizes a perfectly refiecting wall. The 
self-adjointness condition eq. (l4.6p implies 

j(0) = c^(0) W(0) = c^(0)t ° J ^^(0) = c[*i(0)**2(0) + ^2(0)*^i(0)] 
= c[^i(0)*A^i(0) + ^i(0)*A*^i(0)] =0. (4.9) 

In the chiral limit, m = 0, not only the vector current but also the axial current 
is conserved. In the basis we have chosen, the 7-matrices take the form 

1 ^ 1 _ _f 1\ f 1 



(4.10) 
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Thus the axial charge density and the axial current density are given by 

jAix, t) = c^{x, t)"^7°7S^^(a;, t) = c*(x, t)^^(x, t). (4.11) 

Indeed, it is easy to convince oneself that dtPAix,t) + dxjA{x,t) = in the chiral 
limit 771 = 0. Let us now consider the axial current at the boundary x = 

ja{0) = c [|vl/i(0)|2 + |vl/2(0)p] = c (1 + |A|2) |M/i(0)p > 0. (4.12) 

Since in general \E'i(0) ^ 0, the axial current does not vanish at the boundary. 
Hence, chiral symmetry is explicitly broken by the most general perfectly reflecting 
boundary condition. It is well-known that this is indeed the case for the boundary 
condition in the MIT bag model (sj-IEj . Only in the chiral bag model the axial current 
is conserved because it is carried by a pion field outside the bag [l6|. 111 ■ 

It is interesting to ask how the self-adjoint extension parameter A in the rela- 
tivistic case is related to the parameter 7 in the non-relativistic limit, in which 

^2(x) = x^*i(x) = -l-a,^i(x). (4.13) 
2mc zTnct 

Inserting this relation in the relativistic current, 

j{x) = c^(a;)"^a^(x) = c [^i(x)**2(a;) + ^2(x)*^i(a;)] 

= ^[^i(a;)*a,^i(x)-9,^i(x)*^i(a;)], (4.14) 

we indeed recover the non-relativistic current of eq.f l2.2p . In the non-relativistic 
limit, the relativistic self-adjointness condition eg. ( 14. 6 p takes the form 

^ a^.^i(0) = A^i(O) ^ -2mciA^i(0) + (9^^i(0) = 0. (4.15) 



2mci 

Hence, the non-relativistic self-adjoint extension parameter of eg. (12. 46 p can be iden- 
tified as 

7 = -2mczA, (4.16) 
which is indeed real because A is purely imaginary. 



4.2 Reflecting Boundary Conditions for 3-d Dirac Fermions 

Let us now consider Dirac fermions coupled to an external static electromagnetic 
field and confined to a finite domain Q. The corresponding Hamiltonian then takes 
the form 

H = d ■ (j)c + eA{x)^ + f37nc^ — e$(x) = —id ■ Dc + f37nc^ — e$(a;), 
- / (? \ ^ / 1 \ 
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Here 1 and are the 2x2 unit- and zero-matrix, and a is the vector of Pauh 
matrices, while and A{x) are the scalar and vector potential, and e is the 

electric charge. The covariant derivative is given by 



D = y j^i-A{x). 



(4.18) 



The Dirac Hamiltonian acts on a 4-component spinor t). In this case, the 
continuity equation 

5tp(f,t) + V-j(f,t) = 0, (4.19) 

is satisfied by 



p(f, t) = ^(f , ty^^^x, t), j{x, t) = c^'(f, tya'^ix, t). 



(4.20) 



Under time-independent gauge transformations, the gauge fields as well as the Dirac 
spinor transform as 

^$(f) = $(f), '^v4(f) = A{x) - Vifix), ^^(f) = exp ^(^)- (4-21) 

In order to investigate the Hermiticity of the Hamiltonian in the finite spatial 
domain f2, we consider 



d^x x{^) 



d X 



^(f) 

t 



[X) 



-icV + eA{x)j + I3mc^ - e$(x 
-icV + eA{x)^ + (3mc^ - e$(f) x(f) \ \&(x 
— ic dn ■ xi^y (^"^ {x) 

JdQ 

= {^\H\x)*-ici dn-xixyd"^ 
Jan 

which thus leads to the Hermiticity condition 

xixyn{x) ■ a\l/(x) =0, x G dD.. 
We now introduce the self-adjoint extension condition 

A(x) ( ^^5^1 ) , A(f) G GL(2, C), xedn 



(4.22) 



(4.23) 



^3(^) 

^4(x) 

which reduces eq. fl4.23p to 



(4.24) 



xixY 



nix] - a 



n(xi ■ a 



*(f) 



[(Xi(^)*, X2{xy) n{x) ■ aX{x) + {xsixY, XiixY) n{x) ■^]{ ^ 



^i(f) 



'2[X 



X^x) 
Xi{x) 



-nix)-a\ixynix)-a[ ^^fX , , 



(4.25) 
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In order to guarantee self-adjointness of if, i.e. D{H) = D{W), we now demand 

\{x) = —n{x) ■ a\{x)'^n{x) ■ a =^ n{x) ■ crX{x) = — [n{x) ■ crX{x)]'^ . (4.26) 

Hence, n{x)-a\{x) is anti-Hermitean. For Dirac fermions, there is thus a 4-parameter 
family of self-adjoint extensions that characterizes a perfectly reflecting wall. It is 
important to note that the self-adjointness condition eg. (14.241) is gauge covariant 
and implies 



n(x)-j(x) = c^(x)t / _ _^ \m{x) = 

= c[{^i{x)\^2{xY)n{x) ■ a\{x) + (^3(f)*, ^4(f)*) n(f) ■ a] 



^i(f) 



= c{^,{x)\^2{xr) [n{x)-a\{x) + \{x)^n{x)-a] ^ = 0. (4.27) 

Let us again consider the chiral limit m = 0, in which the axial current is also 
conserved. The 7-matrices now take the form 



/ = = ( -1 j • ^' = ^ " = ( j • ^ - ^'' -' '^ V 1 « 

(4.28) 

Hence the axial charge density and the axial current density are given by 

PAix,t) = ^(x,t)V^(^,^), 

l4(f,t) = cvl>(x,t)V77'*(^,t) = -c^ix,ty g ? ) ^(f,t), (4.29) 
such that the axial current at the boundary is 



n[x) ■Ja[x) 



nix] ■ a 



= -c{^,{x)\^2{xr) [n{x) ■ a + X{xMx) ■ aX{x)] |^|;!| ^ .(4.30) 

As in the 1-d case, in general, the axial current does not vanish at the boundary. 
Hence, chiral symmetry is again explicitly broken by the most general perfectly 
reflecting boundary condition. 



4.3 Reflecting Boundary Conditions in the Non-relativistic 
Limit 

In the non-relativistic limit, the lower components of the Dirac spinor are given by 

^■,(x) \ ^' fe + el(f)) / 1 

^- ^ ^ ' ^^^^ ^- a-D^{x), (4.31) 



^4(:?) J 2mc2 V '^2[x) J 2mci 
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with the 2-component Pauh spinor 



(4.32) 



The Dirac equation then reduces to the Pauh equation. Expanding up to the lead- 
ing Zeemann term, but neglecting the higher order spin-orbit coupling and Darwin 
terms, the Hamiltonian entering the Pauli equation takes the form 



H = mc + 



pc + eA{x] 
2mc^ 



e^{x) + na ■ B{x), 



(4.33) 



where B{x) = V x A[x) is the magnetic field and /x = e/2mc is the magnetic moment 
of the fermion. 



Using eq. (14.311) . the relativistic current reduces to 

j{x) = c{^,{xr,^2{xr,^3{xr,^,{xy 



a 



( ^i(^) \ 

^3(f) 
V *4(:?) / 



2mi 
1 

2m 
1 



m{x)^Dm{x) - {Dm{x))^m{x 

-^{x,t)^a X Z3^(f,t) - (Z3^(f,t))1" X a-^{x,t) 
^(f,t)^z3^(x,t) - {Dm{x,t))^^{x,t) 



2mi . 

+ -^V X [^(x,t)^a^(f,t)l . 



(4.34) 



As a curl, the spin term entering the current is automatically divergenceless. In this 
case, the continuity equation 



5tp(f,t) + V-j(f,t) = 0, 
is satisfied with the probability density p{x,t) = \E'(x, t)^\E'(x, t). 
Introducing the gauge covariant boundary condition 



(4.35) 



7(x)\E'(x)+n(^) 
we obtain 
n{x) ■ j{x) = 



D'^ix) -ia X D'^ix) 



0, 7(f) e GL{2, 1 



X e dn, (4.36) 



2mi 
1 

2m 
1 

2mi 



\E'(f)"^n(x) ■ z3\E'(f) - {n{x) ■ D\E'(f))^\E'(f 
^(x, t^fi ■ (^a X D^(x, t)^ -n - (^(Z3^(x, x j ^(f , t) 
[-^(x)S(x)^(f) + ^(f)H(f)^*(x)] = 0, (4.37) 
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which immediately imphes 



(4.38) 



Hence, again there is a 4-parameter family of self-adjoint extensions, now parame- 
terized by a 2 X 2 Hermitean matrix. 

In complete analogy to the previous cases, by partial integration one arrives at 
the Hermiticity condition for the Pauli Hamiltonian of eq. fl4.33p 



dn 



an 



0. 



(4.39) 



One also readily derives 



(jDxi^)^ X ^"^{x) - xi^Ya X Z3^(,f) = V X (x(f)^CT^(x)) . 



(4.40) 



Using Stoke's theorem as well as d{dQ) = 0, (i.e. the boundary of a boundary is an 
empty set), one then obtains 



dfi ■ 



an 



(^Dxix)^ X (T^(x) - x(f)'^CT X Z3^| 



x] 



dn-V X (x(x)^(T^(f)) = / dt- x(x)'f(T^'(f) = 0, (4.41) 
an Jd{dQ) 



such that the Hermiticity condition eq. (l4.39p may be rewritten as 



j dn- (^Dx{x) - ia X Dxix)^ (x) - x{xy (^D"^ { 



x) — iff X D'^ix 



Using the self-adjointness condition eq.f l4.36p . this relation reduces to 



an 



d^x [x(^)W(^)^^(^) - x(^)W(^)^(^)] = 0, 



= 0. 
(4.42) 

(4.43) 



which is indeed satisfied because 7(x) is Hermitean. 

It is again interesting to ask how the Hermitean matrix 7(0;) emerges from the 
matrix \{x) in the non-relativistic limit. Noting that the self-adjointness condition 
eq. fl4.36p can also be expressed as 

7(f)\[^(f) + n{x) ■ aa ■ D'<^{x) = 7(x)\I'(f) + 2mcin{x) ■ a\{x)'^{x) = 0, (4.44) 

one immediately identifies 

7(x) = —2mcin{x) ■ aX{x). (4.45) 
Since n{x) ■ (tX{x) is anti-Hermitean, the resulting matrix 7(x) is indeed Hermitean. 
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5 Domain Wall Fermions 



Using Shamir's variant ^] of Kaplan's domain wall fermions 19|, let us imagine 
that our world has an additional hidden spatial dimension and that we live very 
near a perfectly reflecting fiat domain wall. For simplicity, we first explore this idea 
in (2 + l)-d and then extend it to (4 + l)-d. Since there is no notion of chirality in 
odd space-time dimensions, it is natural to consider massive Dirac fermions. 



5.1 Domain Wall Boundary Conditions in (2 + l)-d 

Let us consider fermions moving freely along the Xi-direction and localized near a 
perfectly reflecting wall located at = at very small positive values of x^. We 
denote the two spatial coordinates by Xi and X3, while time is denoted by Xq. We 
would reserve X2 for Euclidean time which, however, does not play a role in the 
present paper. Starting from the 7- matrices of eq.f l4.10p . it is more convenient to 
change to a chiral basis by performing the unitary transformation 



10/'' ' V 1 



1 \ ::r / 1 \ .^n~. . / 1 



-1 



(5.1) 



The domain wall fermion Hamiltonian then takes the form 

H = ape + f3m(? + S^psc = ape + /3mc^ — ia^cdy,^ (5.2) 
and the 3-component of the current is given by 



c\l/(xi, X3)^a^\l/(xi, X3) 

= ie(*„(.„.3)-,*.(.„.3)-) ( _°i i ) ( ) . (53) 

In complete analogy to the cases discussed before, the Hermiticity condition then 
reads 

l\/^fi(xi,0) 



(x.(x„o)*,x.(x.o)*)^_^ ju:(x;;o) j=o- ^^-^^ 

Introducing the self-adjointness condition 

^«(xi,0) = r/^i(xi,0), (5.5) 
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the Hermiticity condition turns into 



{Xr{xi, 0)* - 0)*r/) ^z.(xi, 0) = ^ xr{xi, 0) = r]*XL{xi. 0). (5.6) 

Thus, in order to ensure D[W) = D[H) and hence self-adjointness, we must demand 

7] = r/* G M. (5.7) 

Since the boundary condition eq. fl5.5p couples left- and right-handed components, 
which transform differently under the (1 + l)-d Lorentz group, it explicitly breaks 
(1 + l)-d Lorentz invariance (unless r/ = or ±oo). 

The Dirac equation takes the form 

pc mc^ + cd3 \ ( -^Rixi^x^) \ _ ( ^i?,(xi,X3 



mc^ - cd^ -pc ) \ ^L(a;i,a;3) ) ^\ ^l(xi,X3) 
Inserting the ansatz 

'^Rixi.x^) = Ajiexp{ipxi)expi->cx3), ^i(a;i,X3) = exp(ipxi) exp(-xx3), 

(5.9) 

for a state localized on the domain wall, the Dirac equation reduces to 



pc mc — cx \ ( Ar \ ( A 



E[ ). (5.10) 



mc + cx —pc J \ Al J \ 
Imposing the boundary condition eq. fl5.5p . one then obtains 



2r] 2 1 — rf I — rf o 2// 

E = — — -mc - — — -pc, cx = — — -mc + — — -pc. (5.11 

1 + 7]^ 1 + 7]'^ 1 + 7]'^ 1 + r/^ 



This is the energy-momentum dispersion relation of a fermion moving with the speed 

I -if 



1 + 7^2 



c < c, (5.12) 



and coupled to a chemical potential 



21] 2 

/i=— — -mc . 5.13 

Only ioT T] = or ±oo, one obtains v = c and /i = 0, as a consequence of (1 + l)-d 
Lorentz invariance. It is important to note that normalizability of the state localized 
on the domain wall requires x > 0, which restricts the allowed range of p. Only 
for ?7 = or ±oo, there is no restriction and one obtains a relativistic massless left- 
handed domain wall fermion with E = —pc. This means that particles (i.e. positive 
energy states) are left-moving [p < 0) while anti-particles are right-moving. 
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5.2 Domain Wall Boundary Condition in (4 + l)-d 

Let us now consider domain wall fermions in (4 + l)-d localized near a perfectly 
reflecting wall located at xs = at very small positive values of X5. We denote the 
four spatial coordinates by xi,X2, X3 and 0:5, reserving X4 for Euclidean time. In this 
case, the 7-matrices are given by eq.f l4.28p . Again, it is useful to change to a chiral 
basis 

(5.14) 

The domain wall fermion Hamiltonian is now given by 

H = a ■ pc + f3m(? + c?p^c = a ■ pc + (3mc^ — ia^cd^, (5.15) 
and the 5-component of the current takes the form 
X5) = c\E'(x, X5)^5^\E'(x, X5) 

- : ) ( ijt::! ) . p^ie) 

where xs) and "^lIxjX^) are 2-component Weyl spinors. The Hermiticity 

condition now reads 

Mf.o)'.x.wo)')(_", :)(i:g:j)-o. 

and the self-adjointness condition takes the form 

^^j(f,0) = r7^L(x,0), r]EGLi2,C). (5.18) 

Inserting this in the Hermiticity condition eq.f l5.17p . one obtains 

{xRix, Oy - xdx, Oyr]) ^l{x, 0) = ^ xnix, 0) = r]^XLix, 0). (5.19) 

In order to ensure D{H'^) = D{H), we must hence demand 

7] = 7]^ = r^ol + ff- a. (5.20) 

In this case, there is a 4-parameter family of self-adjoint extensions. For general 77, 
the boundary condition explicitly breaks (3 + l)-d Lorentz invariance, and for ff ^ 
it even breaks 3-d spatial rotation invariance. 
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The Dirac equation now takes the form 



a • pc mc^ + cd^ 
mc^ — cd^ —a ■ pc 



E 



(5.21) 



In analogy to the (2 + l)-d case, we make the ansatz 



'^r{x,X5) = ARexp{ip-x)exp{-xx5), ^i(f,a;5) = ALCxp{ip ■ x) cxp{-xx5), 

(5.22) 

which reduces the Dirac equation to 



(^3\P\C 

mc^ + cx 



mc 



cx 



- 0-3 IP I c 



A' 



R 



E 



A' 
A' 



R 



(5.23) 



Here we have performed a unitary transformation U{p ) to diagonahze a ■ p, i.e. 



U(p)a . pU{0 = |p |<73, A'^ = U{p)An, = U{p)Al. 



(5.24) 



For simplicity, we now restrict ourselves to the rotation invariant case ff — 0, such 
that the boundary condition reduces to A'j^ — r^o^L with r)o e M. One then obtains 



E 



2^0 2 ^ 
^mc =F 



1 



^01 -I 

p c, cx 



1 "^Imc' ± ^^\p\c. (5.25) 



Again, this is the dispersion relation of a massless fermion moving with the speed 

1 



and coupled to a chemical potential 



1 + %' 

2r;o 
1 + Vo 



c < c, 



T-mc 



(5.26) 



(5.27) 



The normalizability of the domain wall state requires x > which again implies 
restrictions on p. Only for ?7o = or ±00 there is no restriction and one obtains a 
relativistic massless left-handed domain wall fermion with E = \p\c. 



6 Conclusions 



While the results presented here are easy to derive, they seem not to constitute 
common knowledge in quantum mechanics. The theory of self-adjoint extensions is 
not only mathematically elegant, but also physically relevant. Therefore, we hope 
that our paper contributes to changing the view on elementary textbook problems 
such as the particle in a box. While it may not be appropriate to discuss the 
most general reflecting boundary condition in a first encounter with the Schrodinger 
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equation, one might at least point out that other boundary conditions are possible 
as well. 

Since the main purpose of this paper is of conceptual nature, we have not made 
an attempt to determine the value of the self-adjoint extension parameter 7 for ac- 
tual quantum dots. In typical cases, 7 may be very large, so that the wave function 
practically vanishes at the boundary. However, there may be specific situations that 
lead to much smaller values of 7 and thus to quantitatively or even qualitatively 
different behavior, such as bound states localized at the wall of a quantum dot. It 
would certainly be interesting to investigate this in more detail. As a special case 
of the general uncertainty relation for non-Hermitean operators, we have derived a 
generalized uncertainty relation for the self-adjoint position and the non-Hermitean 
momentum operator in a quantum dot. Interestingly, additional boundary terms 
enter this relation. In particular, negative energy states, which may seem to be in- 
consistent with the uncertainty relation, are in perfect agreement with the general- 
ized relation. Minimal uncertainty wave packets, which saturate the corresponding 
inequality and satisfy it as an equality, have been constructed as well. They are 
standard Gaussian wave packets, but require very special boundary conditions and 
thus generically do not exist. If the self-adjoint extension parameter 7(0;) vanishes 
everywhere at the boundary, a constant wave function plays the role of a minimal 
uncertainty wave packet. Furthermore, we have shown that the spectrum depends 
monotonically on the self-adjoint extension parameter 7. 

We have also applied the theory of self-adjoint extensions to theories of relativis- 
tic fermions described by the Dirac equation. In (1-1-1) and (3 -|- 1) dimensions, 
we have considered the most general perfectly reflecting boundary condition as a 
generalization of the one used in the MIT bag model. All these boundary conditions 
necessarily explicitly break chiral symmetry. This can be avoided in the chiral bag 
model were the axial current is carried by a pion field outside the bag. We have 
also discussed generalized domain wall fermion boundary conditions both in (2 -|- 1) 
and in (4 + 1) dimensions. The most general perfectly reflecting boundary condi- 
tion explicitly breaks (1 + l)-d or (3 + l)-d Lorentz invariance or even 3-d spatial 
rotation invariance. A relativistic massless chiral fermion arises only for a particular 
choice of the self- adjoint extension parameters. For other values of the self- adjoint 
extension parameters, one obtains a fermion at non-zero chemical potential, with a 
linear energy-momentum dispersion relation, however, with a speed v < c. 

We conclude this paper by expressing our hope that it may contribute to em- 
phasizing the sometimes subtle differences between Hermiticity and self-adjointness 
also in the teaching of quantum mechanics. 
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